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Abstract 

A hydrodynamic mechanism of interactions of colloidal particles is considered. The 
mechanism is based on the assumption of tiny background flows in the experimental 
cells during measurements by Grier at al. |2|, |3|, |4L pi H. Both trivial (shear flow) 
and non-trivial (force propagation through viscous fluid) effects are taken into account 
for two colloidal particles near a wall bounding the solvent. Expressions for the radial 
(attractive or repulsive) forces and the polar torques are obtained. Quantitative esti- 
mates of the flow needed to produce the observed strength of attractive force are given, 
other necessary conditions are also considered. The conclusion is made: the mechanism 
suggested most likely is not responsible for the attractive interactions observed in the 
experiments of Grier at al.; however, it may be applicable in other experimental real- 
izations and should be kept in mind while conducting colloidal measurements of high 
sensitivity. Several distinctive features of the interactions due to this mechanism are 
identified. 

PACS: 82.70.Dd — Colloids; 47.15. Gf — Low-Reynolds-number (creeping) flows 

1 Introduction 

The motivation for this work comes from the experiments conducted by David Grier at 
al. I], ||, [|, f|, ||, H as well as by other groups [[?], j8|, which show unexplained attraction 
between colloidal particles in suspension. As most other polymer aggregates do, the particles 
used in the experiments dissociate in the solvent with charges of one sign remaining on the 
surface of the particles and charges of the other sign going to the solvent. According to 



the prevailing model for electrostatic interactions of such particles (DLVO || [10], |i~l"f), these 
particles of same-sign charge should repel via a screened Coulomb potential. Several other 
theories for electrostatic interaction mechanisms have been suggested to explain the attraction 
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observed [O, O, [TJ], IT31, but none of them stood the experimental tests. Non-electrostatic 



models have not succeeded in providing concrete reasons for the observed behavior either. 
However, the first hydrodynamic attempt to account for the attraction in the recent work of 
Squires and Brenner |HJ appears quite successful. 

Here we explore another hydrodynamic effect that complements the effect described by 
Squires and Brenner. The nature of attractive force is essentially the same (a non-equilibrium 
hydrodynamic effect based on the method of image forces by Blake Ul7|), but the mechanism 
explored here does not require any charges on the surface of the walls and colloidal parti- 
cles, while the mechanism of Squires and Brenner does. While we do not intend to dismiss 
their mechanism and we agree that it quantitatively captures essential features of the exper- 
imentally observed attraction, we would like to offer a complimentary effect that might be 
responsible for attraction in the cases when the explanation of Squires and Brenner does not 
apply. Examples of such systems are suspensions and meta-stable crystallites, where particles 
are essentially fixed in space (i.e. fluctuate around their equilibrium position). This limits the 
mobility of particles, which decreases applicability of the mechanism of Squires and Brenner, 
but increases applicability of ours (as being "fixed in space" is an essential condition for 
our mechanism). While quantitative description of suspensions and meta-stable crystallites 
requires further work, we feel that the mechanism described below is potentially important 
in accounting for colloidal interactions and complements the picture drawn by Squires and 
Brenner. It also demonstrates a subtlety produced by undetectable hydrodynamic flows that 
is relevant for colloidal measurements. However, this mechanism appears too weak to account 
for the experiments that motivated its development as we currently understand them. 

In this note we first consider physical ideas and assumptions behind this mechanism, and 
then we theoretically investigate two hydrodynamic effects resulting from these assumptions. 
In the discussion section that follows, we estimate orders of magnitude of these effects and 
conclude that the mechanism explored is probably too weak in the particular conditions of 
the experiments of Grier at al. Later we discuss possible other implications and experimental 
tests. 



2 Physical assumptions and geometry 

The experiments of Crocker and Grier were conducted on the samples of colloidal suspension 
confined between a microscope slide and a cover slip. Each experimental run consisted of a 
series of recordings of Brownian motion of two colloidal particles (sufficiently remote from the 
other particles in the suspension) and subsequent analysis of these recordings, which allowed 
one to infer the pair interaction potential Jl], H|. Before each recording the particles were 
caught in the focal plane of the digital video microscopy setup by optical tweezers, which 



are essentially a potential well of the electromagnetic origin jig , , then they were released, 
and their Brownian motion was recorded in several 1/30-of-a-second periods. Subsequently 
the particles were trapped by the tweezers again (in order to be returned to the recording 
area of the sample) and a new cycle began. 

Here we consider the possible effects of background flow in the apparatus. Our motivation 
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originated from the thorough flushes of the solvent carried out with the purpose of removing 
impurities from the solution before each experimental run. However, the flow induced by 
these flushes decays exponentially with characteristic time of the order of d 2 /8i>, where d 
is the distance between the slide and the cover slip and v is the kinematic viscosity of the 
liquid. Hence the flushes cannot be a source of sustained flow: the typical value of this 
time in the experimental conditions is of the order of 10~ 4 -10~ 3 s. At the same time there is 
another effect capable of inducing some flows that are both small enough to be experimentally 
undetectable and large enough to produce the desirable attractive effect during the entire 
experiment (the required magnitude of flow velocity will be shown to be of the order of 
0.1 /im/s). This effect is related to a slight misbalance in pressure in a vacuum system 
connected to the experimental cell. Each cell was connected to the colloid suspension supply 
via two openings in the microscope slide. Colloid refill, ion removal and all other auxiliary 
operations were conducted via this system. When measurements of the attractive potential 
were performed, some negative pressure was applied inside the cell to make the cover slip 
bow inwards, lowing wall separation d from about 50 fivcv to about 5 /iin. While the pressure 
inside the vacuum system was of the order of 15 Torr or 2000 Pa, we estimate^ that the 
difference in pressure between the two parts of the vacuum system had to be only 0.8 Pa in 
order to produce the required magnitude of flow. Thus, 0.04% pressure difference, which is 
definitely beyond the experimental control and can exist unnoticeably in the conditions of 
the experiment [p0|1 , could cause the background flow inside the cell. In general, since the 



required flows are very small, other subtle effects can also contribute to their existence. All 
the subsequent treatment is based on the assumption of existence of small constant flows in 
the experimental cell. One of the purposes of this paper is to show that these undetectable 
flows can lead to quite detectable effects as long as additional conditions outlined below are 
fulfilled. 

In conjunction with this assumption one can be tempted to require that the particles be 
carried along by the flow. Then there would be no hydrodynamic interactions whatsoever 
since particles would not exert any force on the fluid. In contrast, we assume that the 
particles are held in space (by optical tweezers, by inhomogeneities of the wall, by interactions 
with other particles in suspensions or meta-stable crystallites or by some other mechanism) 
while the fluid is flowing by and exerting a Stokes force on them. The mechanism of such 
confinement and its limitations will be discussed in the discussion section. 

Another condition heavily employed below also originates from the experimental pro- 
cedure. All measurements were conducted on colloidal suspensions enclosed between two 
microscope slides separated by a small gap (of the width of a few microns) so that the par- 
ticles were confined in the vicinity of at least one of the slides. Moreover, the strongest 
attraction was observed when particles were located sufficiently close to one of the slides. 
Therefore, we consider particle interactions near a plane wall. As it will be shown, both 



lr rhe only non-trivial component of the Navier-Stokes equation reads ^ = ?7§jt, with x being along the 
flow and z being normal to the walls (77 is dynamic viscosity). Since pressure p does not depend on z, the 
standard solution is the parabolic profile for velocity u — r yz(d — z)/d, with d being the distance between the 
parallel walls. This yields \Ap\ — 2rij\Ax\/d for pressure difference between the openings separated by Ax. 
Typical values Ax = 2 cm, 77 = 1 mPa ■ s, and d = 5 fim give a result |Ap| = 0.8 Pa for 7 = 0.1 s — 1 . 
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Figure 1: Geometry of the problem. 



effects described below are absent in the infinite space. 

The last assumption used will be the absence of inertial effects or their negligible con- 
tribution, which is a reasonable approximation for the case of small Reynolds numbers (the 
numbers involved are of the order of 10~ 5 -10~ 6 ). Gravity is also unimportant due to roughly 
equal densities of the material of the particles and the fluid, and so the only forces acting 
on the particles are of the hydrodynamic origin. Thus, we consider the hydrodynamic inter- 
actions of spatially fixed particles subject to slow flows (with small Reynolds numbers) in the 
semi-infinite space. 

In particular, consider particles 1 and 2 in the flow field above a plane wall. Let particle 
1 be at height hi above the plane and particle 2 be at h 2 , and let the full (three-dimensional) 
distance between the particles be b (Fig. [l]). In agreement with the experimental conditions, 
we picture both particles as spheres; let their radius be a. We do not assume any relation 
between h and b, in particular they can be comparable or one of them can be much larger than 
the other. However, we do assume that the size of the particles a is much smaller than any 
other distance scale (h or b), so that the particles can be thought of as point-like. This last 
assumption is not strictly obeyed in the experiments (although particles are smaller than the 
separations between the objects), but it greatly simplifies the consideration and presumably 
does not change the qualitative character of the interaction potential. 

We choose the origin of the system of coordinates at the location of particle 1. Axis 
z is directed away from the wall, so that the wall is defined by the equation z = —hi. 
Axis x is directed along the velocity of the flow u (parallel to the wall). Particle 2 has 
spherical coordinates (b,9,<f>) or Cartesian coordinates (b sin 9 cos 0, b sin 9 sin 0, b cos 9) with 
the standard choice of angles 9 and 0. Then hi = hi + bcos9. (See Fig. [1] for geometry.) 



3 The shear flow effect 

The simplest effect in the system described above comes from the flow gradient near a plane 
wall (see e.g. Landau [21]). Sufficiently close to the wall the flow speed u must grow linearly 
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with the distance to the wall: u = 7/ix, where 7 is the shear modulus. Then the original 
flow field at the location of particle 1 creates a Stokes force on that particle Fi = 6717707/110;, 
while at the location of particle 2 this force is F2 = 6717707/122:, where t] is the fluid viscosity. 
Obviously, the higher particle experiences a higher force from the fluid flowing by. Now, if we 
define the apparent interaction force between the particles as the difference between the radial 
components of forces acting on each particle: F r = F 2 • f — F x • f (so that the positive force 
corresponds to repulsion and the negative one corresponds to attraction), then the difference 
in forces exerted on the particles can be interpreted as either repulsion (if particle 2 is in the 
1st or the 3rd quadrant of the xz-plane with respect to particle 1) or attraction (if particle 2 
is in the 2nd or the 4th quadrant). The exact result for the defined above interaction force is 

F r = 6^77706 cos 4> cos 9 sin 9. (1) 

Of course, the interpretation of this force as repulsion or attraction does not mean that one of 
the particles acts on the other, but in the conditions of the experiment (where the potential 
is inferred from the measurements of Brownian motion for a fraction of a second) there is 
no simple way to distinguish whether particles interact directly or just move under the two 
external forces as if they interact. 

If the two particles are at the same height above the wall, then no apparent interaction 
force is present. This can be seen both from expression ([!]) and from the fact that the Stokes 
forces exerted on both particles are exactly the same. It can also be shown that there is no 
torque or force making particles leave the original height. 

Expression ([[]) allows one to estimate the minimal value of the flow velocity necessary 
to produce the observed magnitude of the interaction force. The maximum of the absolute 
value of angular dependence of the force is reached at 9 = 7r/4 or 37r/4 and <fi = or 7r; 
at this configuration | cos^cos^sin^l = 1/2. Taking the typical values F r = 1CT 15 N, a = 
0.5 x 1CT 6 m, b = 2 x 1CT 6 m, 77 = 1CT 3 Pa • s |], ||, |2(|, one can easily obtain: 



7>0.1s~ 1 . (2) 

Thus, flows as low as 0.1 /im/s (at height of 1 fim above the wall) can produce the ob- 
served magnitude of interactions. These values of the velocity are below the experimentally 
detectable level, and such flows can exist in the system [pOfl . 

Now, if there are two alternatives (attraction in 2nd and 4th quadrants and repulsion in 
1st and 3rd quadrants), why would only attraction be observed? A possible explanation will 
be deferred until the next section. 



4 Force propagation via viscous fluid 

Apart from the trivial effect described above, there is a more elaborate mechanism of particle 
interactions via viscous force propagating through the solvent. If a point force F ( "stokeslet" ) 
is applied at point x, then the velocity field perturbation v at point y is linearly related to 
the stokeslet magnitude: fj(y) = Hijiy — x)-Fj(x), where Hij is known as the Oseen tensor 
(in the case of an extended source of force an integration over coordinates is required, so 
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that the Oseen tensor is just a Green function for the velocity field). The force exerted on 
a spherical particle at point y is then a Stokes force created by the velocity perturbation v: 
f = 6irar]v. In the case of two particles in viscous fluid the Stokes force Fi exerted by the 
original flow on particle 1 (given in the previous section, but taken with a minus sign, since 
now the particle, being held in place, acts on the fluid) will be the source of the perturbation 
of the velocity field v 2 at the location of particle 2 and hence the source of the perturbation 
force f 2 acting on particle 2, and vice versa. Thus, the interaction force between the particles 
can again be defined as the difference between the radial components of forces acting on each 
particle: f r = f 2 -f — fi - f, but now this will be the true interaction force, the one propagating 
through the fluid. 

Expression for the Oseen tensor Hij in the infinite space is well known (see e.g. Happel 
and Brenner ||22|| ) and leads to the identical zero for the radial interaction force (the flow is 
uniform). Thus, no hydrodynamic interactions are possible in the infinite space. The Green 



function Hy for the semi-infinite space was constructed by Blake [17|, although it can be 



alternatively derived by the method of Lorentz ]2"3|]. Calculation of the interaction force f r 

(3) 



based on his result leads to the following expression: 

1 + U x t 2 + 6tit 2 (ti + t 2 ) 



f r = 97r?77a cos <b cos 9 sin 6 ... 

(1 + 4M 2 ) 5 /2 

where t± = h±/b and t 2 = h 2 /b = t± + cos8. The asymptotics of the above result are: 

h 

f r = 547r7?7a 2 — cos cos 2 9 sin 3 9 if hi <C b (4) 

and 

f r = 97r?77a 2 cos <p cos 9 sin 9 1 — — J if h 3> b, (5) 



where in the last line h can be either hi or h%. Typical behavior of f r as a function of 
l/ti —b/hi for 9 = it/A and <fi = is shown on Fig. |^. 

Inspection of the result (|3|) indicates that the coefficient preceding this force is a 2 , in 
contrast to the factor ab in eq. (|]). However, far from the wall the angular dependences of 
the kinetic force ([I]) and the real force ([$]) are identical. Thus, no new effect is present here: 
although the force magnitude is different, the signage is the same — attraction in 2nd and 
4th quadrants and repulsion in 1st and 3rd quadrants. So, expression (£J) serves just as a 
correction to the main result ([[]) (since, by assumption, a <C 6). It also gives the same order 
of magnitude for 7 as estimate (0) does. 

Note also that f r of eq. @ reinforces F r of eq. ([!]). This is easy to understand from 
the fact that the perturbation velocity v (and hence the perturbation force f) is in general 
oppositely directed with respect to the original Stokes force F acting from the fluid on the 
particle, or similarly directed with respect to the force — F from the particle on the fluid (see 
Fig. ^ for the infinite space), and is linearly proportional to the magnitude of the force F. 
Thus, for instance, if the radial projection Fi • r is greater than F 2 • f, then the negative 
radial projection — f 2 • f should in general be greater than — fi • f , leading to the same signage 
of the above defined F r and f r . 
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Figure 2: Interaction force f r /{rj^a 2 ) as a function of distance bjh\ for = 7r/4 and = 0. 



-F 



u 



Figure 3: Typical orientation of bulk flow velocity u, Stokes force exerted by the fluid on the 
particle F, Stokes force exerted by the particle on the fluid — F, and velocity perturbation 
field v(x) at some distance b from the particle (infinite space). 
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While the direct effect of the true hydro dynamic forces is relatively weak comparing to 
the kinetic effects of the previous section, they have an indirect effect that may explain the 
tendency towards observable attraction. In particular, if the two particles are at the same 
height above the wall, then there is a force in vertical direction pushing each of the particles 
into the 2nd or the 4th quadrant with respect to the other particle: 

t 4 

f z = -547T77 7 a 2 cos0 (i + 4t2)5/2 , (6) 

where t = h/b and h = hi = h 2 . For h = b, <fi = 0, 7 = 7 m i n = 0.1 s -1 and the same typical 
values of a and 77 as in estimate (|2|) the absolute value of this force is 7.6 x 10~ 17 N, which 
is about 8% of the typical value of F r . In general, for arbitrary h\ and /i 2 , this effect can be 
described by the torque acting in the ^-direction, which does not vanish at 9 = n/2 (unlike 
in the effect of the previous section). This torque is defined as T = (f 2 • 9 — fi • 9)b/2, and a 
general expression for T can be obtained on the same grounds as result (R|): 



—rj'ya ocos0 



, v l + 4M2) cos 2 9 + 6* 1 i 2 (i 1 + t 2 f cos Id 



cos' 9 - 



;i + 4m 2 )5/2 



(7) 



Thus, this torque indeed favors the attraction geometries by pushing particles into the config- 
urations where attraction is created by the shear flow effect of the previous section. Therefore, 
attraction indeed should be observed more often than repulsion if particles are positioned at 
the same height above the wall at the beginning of each experimental run. 



5 Discussion 

While the mechanism described in the previous section seems to lead to the correct results 
and be a feasible explanation to the observed attraction, there remains one question: What 
holds the particles fixed in space? As we mentioned earlier, the difference in fluid and particle 
velocities is the necessary condition for the existence of the Stokes forces producing all the 
above effects. On the other hand, even if a particle is initially set at rest in the flow, while 
not being held by some external force, it acquires the velocity of the flow for relaxation time 
r (defined by v (t) = u(l — exp(— t/r))) much shorter than the duration of each recording 
period (1/30 s), and hence it does not exert any force on the fluid for most of the recording 
time. For instance, for a spherical particle of radius a = 0.5 x 10~ 6 m and of density of the 
water (p = 10 3 kg/m 3 ) in a fluid of viscosity of the water (rj = 10~ 3 Pa • s) this relaxation 
time is 

r= 2P^ 56ns _ (8) 
9 77 

Thus, the described effects require an external force that holds the particles at a non-zero 
velocity with respect to the flow. 

The most probable candidate for this role could be a pair of optical tweezers used to trap 
the particles before each recording run. One scenario is straightforward to consider. While 
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being held by the tweezers before each recording of their Brownian motion, particles are 
positioned to an attractive configuration by means of the torque described in the previous 
section. After the holding potential of the tweezers is removed and the recording starts, the 
dominant shear flow effect influences the Brownian motion of the particles, and an attractive 
"potential" is registered. Note that this does not even require any Stokes forces to be present 

— just a drift in the flow with different velocities (so that the further from the wall particle 
moves faster) influences the Brownian motion in essentially the same manner. Since the 
initial alignment was practically always chosen in the focal plane (so that both particles were 
at the same height above the wall) and along the longest dimension of the slides (so that 
one of the particles was always down the flow with respect to the other) [^0| , one might 
speculate that the hydrodynamic effects proposed above had to produce apparent attraction 
in practically all the conducted experiments. 

However, in spite of general feasibility of such a mechanism, its estimated magnitude is 
discouraging. In the particular conditions of the experiment the force gradient in the optical 
tweezers is of the order of 1 pN//im = N/m pCfl , which leads to the forces of 3 or 4 
orders of magnitude higher than force (J6|). Thus, particles are not easily displaceable by the 
hydrodynamic forces while being held by optical tweezers, and therefore the hydrodynamic 
mechanism appears too weak if the tweezers work as it is commonly understood. (Note, 
however, that the estimate for the force @ was made for the minimal predicted value of the 
shear modulus; higher flows lead to stronger hydrodynamic interaction.) Therefore, optical 
tweezers, in spite of being a promising candidate, most likely cannot perform the role of the 
required external force. 

Other candidates to this role might also be possible (e.g. forces due to the charges on 
the wall, due to inhomogeneities of the wall, or concerted mean-field-like forces at high- 
concentration systems like suspensions and meta-stable crystallites); however, we were not 
able to find a feasible scenario of how exactly these forces would keep the particles fixed (or 
at least at non-zero relative velocity with respect to the flow) . Thus, while the hydrodynamic 
mechanism certainly exists when such forces are present, the main difficulty in possible ac- 
counting for the attraction observed in the experiments like those of Grier at al. comes from 
identifying these forces. At the same time, in the situations where there are such forces, the 
hydrodynamic mechanism considered here is always at work (and this should be kept in mind 
while conducting high-sensitivity force measurements in colloidal systems). One example 
where the required forces exist can be the systems of meta-stable crystallites [||, ^3). A self- 
sustaining mechanism can be employed here: The attractive force holds particles together, 
creating thereby the necessary condition for the described effect (flow past fixed particles 

- assuming some flow is present), which in its turn leads to the existence of attraction. 
Further work would be required to account in detail for such a self-sustaining mechanism. 
Here we only notice that while the concerted forces acting on a particular particle can be 
mean-field-like, the whole system of particles cannot be held at rest by only the interaction 
forces between the particles in that system, as in this case the whole system would be carried 
along by the flow. Therefore, there should be some external forces different from interactions 
with the neighboring particles and acting at least on the boundary layer of the crystallite 
system. 
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As it is also apparent from our consideration, flow is not the only possible source of the 
described effect. What is necessary is the force exerted by the particles onto the fluid. This 
can be achieved not only by making a fluid flow by the fixed particles but also by making 
particles move in the stationary fluid. Thus, for instance, particles can be dragged by some 
external force in such a way that the further from the wall particle experiences higher force; 
this should lead to effectively the same result for the interaction potential. Essentially, this 
has a lot in common with the mechanism of Squires and Brenner where particles are 
dragged away from the wall by the electromagnetic repulsive force. 

Note that presence of a second wall, positioned symmetrically with respect to the interact- 
ing particles, leads to canceling of all considered effects (both force and torque). Therefore, 
another restriction on the possible experimental geometries where this hydrodynamic mecha- 
nism could be responsible for interactions is that the particles must be located away from the 
center plane of the suspension sample confined by the microscope slides. In practice, however, 
this restriction is easily avoided due to a slight difference in particle and fluid densities, so 
that the particles are pulled closer to the lower wall by the gravity (if their density is higher) 
or pushed closer to the upper wall by the Archimedes force (if their density is lower). 

The proposed mechanism immediately suggests several possible experimental tests on the 
presented effects of flow. Most of these tests originate from the following list of the predicted 
properties: 

1. The simplest hydrodynamic effect originates from the velocity gradient near a surface 
and can lead to the apparent attraction or repulsion (although no direct interaction 
exists). Both outcomes are equally probable with randomly chosen initial geometries. 

2. The true hydrodynamic interactions propagating through the viscous fluid do exist, 
but they lead to the same qualitative dependence on the geometry — attraction in 2nd 
and 4th quadrants and repulsion in 1st and 3rd quadrants. Nevertheless, there exists 
a mechanism favoring attraction geometries. 

3. The flows necessary to produce any hydrodynamic effects can be as low as 0.1 /xm/s 
(more generally, 7 > 0.1 s _1 ). 

Different initial geometries fixed by optical tweezers can be tried to check if the interaction 
strength changes depending on the position of the particles. This requires only relocation 
of principal potential minima of the tweezers, which can be relatively easily achieved by 
refocusing or rotating them. If both attractive and repulsive results are recorded (or if a 
stronger interaction produced by some other mechanism gets modified accordingly), then the 
effect is at work. 

One can also look explicitly for the flow. This flow should be easily detectable if particles 
are allowed to drift freely for a sufficient interval of time. The flow could also be deliberately 
induced at higher levels to increase the magnitude of the hydrodynamic force and check that 
the mechanism works as predicted. 

The mechanism discussed here and the method to estimate its magnitude were proposed by 
T.A. Witten, who also provided ongoing advice during this work. D.G. Grier as well as other 
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people of his group (especially J. Plewa and S. Behrens) supplied essential experimental input and 
a number of ideas about the proposed mechanism and its feasibility. This work was supported in part 
by the MRSEC Program of the National Science Foundation under Award Number DMR-9808595. 
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